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Abstract 

We discuss general 2-fluid hydrodynamic equations for complex fluids, where one kind is 
a simple Newtonian fluid, while the other is either liquid-crystalline or polymeric/elastomeric, 
thus being applicable to lyotropic liquid crystals, polymer solutions, and swollen elastomers. 
The procedure can easily be generalized to other complex fluid solutions. Special emphasis is 
laid on such nonlinearities that originate from the 2-fluid description, like the transport part of 
the total time derivatives. It is shown that the proper velocities, with which the hydrodynamic 
quantities are convected, cannot be chosen at will, since there are subtle relations among them. 
Within allowed combinations the convective velocities are generally material dependent. The 
so-called stress division problem, i.e. how the nematic or elastic stresses are distributed between 
the two fluids, is shown to depend partially on the choice of the convected velocities, but is 
otherwise also material dependent. A set of reasonably simplified equations is given as well as a 
linearized version of an effective concentration dynamics that may be used for comparison with 
experiments. 
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1 Introduction 



The thermodynamic and hydrodynamic properties of multi-component complex fluids are 
determined by the microscopic degrees of freedom of their constituents and the coupling 
between these degrees of freedom. Such systems can exhibit rather rich phase behavior 
and dynamics, especially when one or more components is a structured or macromolecular 
fluid Pp. Due in part to the coupling of internal degrees of freedom, these systems can 
also exhibit novel flow-induced structural evolution phenomena, including shear-induced 
phase transformations and flow alignment of constituents on microscopic to mesoscopic 
length scales. Such structural evolution in turn leads to nonlinear rheological behavior, 
such as stress overshoots in response to imposed rates of strain, plasticity, and thixotropy. 

The overwhelming complexity of the microscopic description of these systems, such 
a detailed description is often not well suited for analysis of the macroscopic dynamical 
behavior. Instead, explicit macroscopic models have been developed for this purpose. 
Some such models have been obtained by a suitable coarse-graining procedure starting 
from a microscopic theory. Others are purely phenomenological models constrained only 
by conservation laws, symmetry considerations and thermodynamics. The so-called "two- 
fluid" models for binary systems of distinct components or phases are useful examples of 
such a macroscopic approach [2]. In the two- fluid description, each component or phase 
is treated as a continuum described by local thermodynamic variables (e.g. temperature, 
density, and relevant order parameters), and dynamical quantities (e.g. velocity or mo- 
mentum). In general, these variables for the constituents are coupled. For instance, the 
effective friction between components in a binary fluid mixture leads to a drag force in 
the macroscopic description that is proportional to the local velocity difference. 

Two-fluid models have been employed in many different physical contexts. The two- 
fluid approach is a key element of many traditional models for multi-phase flow of bubbly 
liquids, fluid suspensions of particulates, and binary mixtures of simple fluids [3]. Other 
examples in condensed matter physics include two- fluid models for superfluid helium 
dynamics of plasmas [5 , transport in superconductors [Uj, viscoelasticity of concentrated 
fluid emulsions [7], flow-induced ordering of wormlike micelle solutions [S], flow of colloidal 
suspensions [9, . Two-fluid models have been used extensively to model a wide range of dy- 
namical phenomena in polymer solutions and binary blends, including the hydrodynamics 
modes of quiescent polymer solutions [HUE], kinetics of polymer dissolution [T2|, hydro- 
dynamics and rheology of polymer solutions and blends PBJ-JIII, an d polymer migration 
and phase separation under flow [2T)]-j2*T]. 

These examples share certain general features. In each, two distinct species or coex- 
isting phases (gas and liquid, normal fluid and superfluid, polymer and solvent, mesogens 
and solvent etc.) with mass densities p\ and P2, which are conserved individually in the 
absence of chemical reactions, move with distinct velocities V\ and t>2, respectively. Due 
to (usually strong) internal friction, the momenta of the constituent species, p\V\ and 
P2V2, are not conserved individually. Of course, total momentum is conserved. In most 
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cases of fluid mixtures the friction is so strong that the velocity difference Vi — v 2 is 
nonzero for very short times only, i.e. it is a very rapidly relaxing quantity that is not 
included in the hydrodynamic description for binary mixtures. However, there are sys- 
tems and situations, where the relaxation of the relative momenta is slow enough to have 
a significant influence even on the hydrodynamic time scale. Then a two-fluid description 
is appropriate and useful. 

In this communication we focus on a general nonlinear two-fluid description of complex 
fluids, where one species is a viscous Newtonian fluid and the other either a polymer or 
a liquid crystal. Emphasis is placed on the rigorous derivation of the dynamic equations 
within the framework of hydrodynamics as contrasted to ad-hoc treatments. The resulting 
equations are rather general and complicated. They can and have to be simplified for 
special applications or systems by appropriate and well-defined approximations. One 
of the advantages of starting from the general theory is the possibility to identify and 
characterize the approximations made. The hydrodynamic method, described in some 
detail in j^Hj-jSO], is quite general and rigorous, being based on symmetries, conservation 
laws, and thermodynamics. In the following sections, we provide a detailed analysis of 
two-fluid models for lyotropic nematogens in a simple viscous solvent, followed by an 
abbreviated extension of this treatment for isotropic elastomers (e.g. entangled polymer 
solutions and gels) in a simple viscous solvent. We close with a discussion of our general 
results and their possible implications for experiments. 



2 Thermodynamics 

The hydrodynamics of fluid mixtures as described above is governed by conservation 
laws (individual masses, total momentum and total energy), balance equations for the 
liquid crystalline degrees of freedom, for the transient elasticity of polymers and for the 
relaxation of relative momentum. There are different ways of writing the appropriate 
equations. One popular choice is to use equations for individual mass densities and 
individual momentum densities, another to use the mass density and one concentration 
variable and the total momentum density and the relative velocity difference. Since they 
both have their advantages and disadvantages we will present both ways of description 
and show, how they are connected. In this and the following sections we will use a nematic 
liquid crystal as the second, complex fluid. Transcription of the formulas to the polymer 
case will be given in Sec.flHJ). 

The starting point of any macroscopic description is the total energy E of the system 
as a function of all the relevant variables. Since the energy is a first order Eulerian form 
of the extensive quantities, we can write 

E = eV = J edV = E(M U M 2 , V, G x , G 2 , S, M 2 V M 2 5n i ) (1) 

The masses, Mi, M 2 and momenta G ly G 2 of species 1 and 2 are related to the appropriate 
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(volume) densities by p 1 = M 1 /V, p 2 = M 2 /V, g 1 = G x /V = p 1 v 1 , g 2 = G 2 /V = p 2 v 2 , 
while for the entropy density a = S/V. The nematic degrees of freedom are related to 
species 2 and consist of director rotations 612^. The nematic phase shows orientational 
order along the line denoted by n (with n 2 = 1) called the director. Since up and 
down (along that line) cannot be discriminated, all equations have to be invariant under 
a n — > — n transformation. Homogeneous rotations do not cost energy, so in a linear 
description (of the field- free case) 5rii is absent in E and only gradients enter j3*T| . 
We have kept both terms to cope with the general case. 

Introducing thermodynamic derivatives (partial derivatives where all other variables 
are kept fixed) we define temperature T, thermodynamic pressure p, chemical potentials 
Hi, p 2 and velocities v x , v 2 of the two fluids, as well as the conjugate fields Xij an d ki 
connected to the nematic degrees of freedom 

dE de dE de dE de 

Mi = ^TTT = 7\ j fJ-2 



dS da ^ dMt dpi ^ A dM 2 dp 2 

dE dE de dE de 

P= ~dV> "^dCr^ V2= dG 2 = dg- 2 (2) 
dE de dE de 

* ~ d(M 2 m) ~ d{p 2 m) ' Xij ~ d{M 2 V jni ) ~ d(p 2 V 3 ni) 

Expanding eq. into first order differentials, the condition dV = leads to an expression 
for the pressure 

p = -e + Tcr + p^ + p 2 fx 2 + v 1 ■ g 1 + v 2 ■ g 2 (3) 

where we have introduced the effective chemical potential of the nematic fi 2 = pi 2 + 
Xij^jUi + kibrii. In addition, the differentials are related by the Gibbs relation 

de = Tdo + fii dpi + p, 2 dp 2 + v% ■ dg 1 + v 2 ■ dg 2 + *&ij dVjTii + hi drii (4) 

with the more familiar nematic conjugate fields = p 2 Xij an d hi = p 2 ki. From eqs.(JHl 
HJ) the expression for the differential pressure results (Gibbs-Duhem relation) that is useful 
in switching from pressure to chemical potentials or vice versa 

dp = a dT + pi dpi + p 2 dp 2 + g x ■ dvi + g 2 ■ dv 2 — dVjrii — hi drii (5) 

A second set of equations is obtained by switching to the total density, p = Pi + p 2 , 
and the total momentum, g = g x + g 2 = p x v x + p 2 v 2l which are the sums of the original 
quantities and which are both conserved quantities. The two-fluid nature has then to be 
represented by additional variables. A natural choice seems to be the use of the density 
and momentum differences. However the latter choice is problematic, since it necessarily 
implies the conjugate quantities also to be the (arithmetic) sums and differences of the 
original conjugate quantities. Thus, the conjugate to g would be v x + v 2 , which does 
not reflect correctly the possible one- fluid limits p x — > or p 2 — > 0. The physically 
acceptable conjugate to the total momentum is the mean velocity v defined by p~ x g. 
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Insisting on v, the mean velocity, to be the conjugate of the total momentum g, the 
choice of the remaining variable describing the different velocities is severely limited. 
Compatibility with (j3J) allows as variable only the velocity difference *** w = v± — v 2 
(with m = p~ x p\p 2 w as conjugate quantity) or more generally aw as variable with 
a~ l p\p 2 p~ l w as conjugate, where a can be freely choosen. There is no a-priori advantage 
for any of the choices and we will stick to a = l. 1 *** From w = g 1 / pi — g 2 / p 2 one gets 

v\ = p^g + (1 - <p)w, v 2 = p~ 1 g-(f)w (6) 

The representation of the two different densities is less problematic. A convenient choice 
for that variable is the concentration, = pi/p, with p 2 /p = 1 — 0. If the expansion 
coefficients of the two fluids are the same, can be interpreted as the volume fraction as 
well. Instead of one could have used, e.g. the density difference p\ — p 2 (or any other 
linear combination of p\ and p 2 different from p) as variable without much change. 



After some trivial algebra eqs.(JB]EI) can be written in the new variables as 

p = -e + Ta + pn + p^g 2 (7) 

de = T da + U' d(f) + fx dp + v ■ dg + m ■ dw + dVjUi + hi drii (8) 

dp = a dT + p dfi + g ■ dv — m ■ dw — H'd(f)— dVjUi — hi drii (9) 



where we have introduced the relative pressure II', the total chemical potential /i, the 
mean velocity v and the weighted relative momentum m defined by 

IT = p(fi 1 -fi 2 ) + w-g + pw 2 (l-2(j)) = pli 

H = /ii0 + /i 2 (l - 0) + W 2 (j)(l - 0) 

or vice versa p,\ = p, + p _1 p 2 (II — w ■ Vi) (10) 

fx 2 = p, — p" 1 p\ (11 + w ■ v 2 ) 

where v = <pv 1 + (1 - 0) v 2 = p _1 (c/i + g 2 ) 

m = p(l -(f>)(f>w = (p 2 g 1 - pig 2 )p~ l 

The Gibbs relations connects variables that show different rotational behavior. Energy, 
entropy, the densities and the concentration are scalar quantities that do not change un- 
der (rigid) rotations, i.e. de = da = dp = dp x = dp 2 — dtp — 0. The vectors and 
tensors are transformed according to drii = Qijfij, dgi = Qijgj, dwi = QijWj, dVjUi = 
^jfcVfcnj + VjQikTik, where Qij = —flji is any constant antisymmetric matrix. The rota- 
tional invariance of the Gibbs relation (J4I8|) then leads to the relation 

hiUj + ^kiS/jUk + ^>i k V k n,j = hjrii + ^ kj Vin k + ^ jk V k ni (11) 

which has to be fulfilled by the conjugate quantities. There are no contributions from the 
momenta and velocities, since g || v, w \\ m, and g lj2 \\ v 12 . Relation (fTTj) is useful for 
reformulating the stress tensor, in particular to symmetrize it explicitly 

1 The choice a = P1P2P 1 would just interchange the roles of w and m as variable and conjugate. 



5 



Having set up the thermodynamics of the relevant variables we are now in a position 
to establish the structure of the dynamic equations. 



3 Dynamic Equations 

For the two fluids there are independent continuity equations stating that neither mass 
can be destroyed nor created, but only transported. Transport can involve convection as 
well as (relative) diffusion. This leads immediately to 

Pi + Vi(piv i m +j®) = (12) 
P2 + V l (p 2 v? ) -jV) = (13) 

When dealing with components of vectors, the subscripts 1, 2 are written as superscripts 
for clarity. The phenomenological mass currents in Eqs. (jl2l [T3*j) add up to zero, since the 
total mass current is equal to the total momentum density g = pv. Eqs. (fT2l fT3*j) can be 
rewritten in terms of the total density and the concentration as 

p + Vjpvj = (14) 
+ ^V,0 + p- 1 V,(p0(l-0M+j/ 1) ) = (15) 

show the characteristic difference between extensive quantities, where convection is of the 
form V ■ (v*) and intensive ones with v ■ V*. 

Note that the concentration does not obey a conservation law, except when linearized 
around a zero- velocity state or if p = const, is assumed. Because the mass current density 
of the total fluid is equal to the momentum density g (= pv), the total mass is convected 
by the mean velocity in (|14j) . In Eqs. (|T2| HHJ) the convective terms are not fixed 
a priori, since the phenomenological current j[ can contain contributions proportional 
to some velocities, thus altering the effective velocity, with which the different quantities 
are convected. We will discuss this point extensively after having derived the full set of 
equations. 

The dynamic equations for the other variables are even more complicated and also 
contain phenomenological parts. These are expressed by yet to be determined currents. 
But they also contain convective (or transport) terms. Therefore, we can set up the 
following equations as an ansatz 

(16) 
R/T (17) 
(18) 
(19) 
(20) 



e + V j (e + p)v j + V lJl (e) = 

& + Vjtrvj + V i j i (f7) = 

9% + Vj9i v j + Vj-Oij = 

■ibi + VjVjWi + Xi = 

hi + VjV jfii + Yi = 
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containing either the divergence of a current (j i , j^, j/ , a^) when conservation laws 
are involved, or quasi-currents (Xj, Yj) in the case of balance equations for non-conserved 
variables. Each of the currents and quasi-currents consists generally of three parts: A 
geometric or symmetry related one without any phenomenological coefficients, which we 
will determine below, and two phenomenological parts, which are either reversible (super- 
script rev) or irreversible (superscript dis). The phenomenological parts will be discussed 
in section |21 The entropy balance (|17|) is not a conservation law, since for irreversible pro- 
cesses the entropy production R has to be positive and only for purely reversible actions 
R = 0. In Eqs. (jl4H 20 j) the convective terms are written down such that all quantities 
are convected by the same velocity 2 This is dictated by the postulation of zero entropy 
production (these transport terms are reversible). However, it should be repeated that 
the phenomenological reversible currents may change the effective convection velocity, 
something we will discuss later. 

Putting the dynamic equations (JT3HHI1) into the Gibbs relation (JU) the condition R = 
(R > 0) for the convective and the reversible (dissipative) phenomenological parts of the 
currents, leads to the following conditions 





= 5 ijP + V 4 n fc + + of s) 


(21) 


x t 


= ViU + x} rev) +x l (dis) 


(22) 


Yi 
(*) 

i 

;« 


_ y ( rev ) _j_ y ( dis ^ 

■ (a,rev) , ■ (aAis) 
= Ji +Ji 

■ (l.rev) , ■ (l,dis) 
= Ji +Ji 


(23) 
(24) 
(25) 



with the generalized conjugate to the nematic degrees of freedom hi = hi — Vj^y = 
P2{ki — V jXij)- The stress tensor <7y contains the isotropic pressure p (J7J), while the quasi- 
current Xi of the relative velocity contains the gradient of IT, the relative pressure divided 
by the total density, (|10p. The terms related to the nematic degrees of freedom are well- 
known from ordinary nematodynamics. The energy conservation law is redundant here, 
because of the Gibbs relation (j3J) and is not needed. 

The phenomenological parts have to fulfill (up to an irrelevant divergence term) 3 

R = -j^ViT - j^ViU - agVjVi + fu r> } + m t X t { * ] > (26) 

with the equal sign (> sign) for * = rev (* = dis), respectively. 
Eq. (j26|) also reveals the equilibrium conditions 

ViT = V 4 n = Aij = hi = mi = (27) 

where 2Aij = VjVi + ViVj. 

2 In the energy conservation law l|16l) the free enthalpy e + p is convected, cf. 30 
3 The true condition is / R dV > 0. 
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Before we will determine the phenomenological parts in (|21H24j) . we first have a look 
into the 2-fluid statics. 



4 Statics 

The statics is given by the connection of thermodynamic conjugates with the variables. 
The conjugates are defined by partial derivatives of the energy density 0- Thus one 
can either write down a phenomenological energy expression and take the derivatives or 
give directly these relations under the proviso that mixed derivatives are equal. Two of 
these connections have already been given in eq.()10|) relating v with g and m with w. Of 
course, these are not really static relations. They are fixed (and not of phenomenological 
nature), since the mass current pv is identical to the momentum density and since the 
kinetic energy density is {l/2)p\v\ + (1/2)^2^2- 

The 3 scalar conjugates {T, II, p} have to be expressed by the variables {a, 0, p} or 
using the other set of variables {a, p 1 , p 2 } and conjugates {T, p 1 , p 2 } by 

5T = ^ r 5a+—5p 1 + —5p 2 (28) 

Cy pOli P&2 

Pi = -^—5p 1 + ~^—5p 2 + —5a (29) 
p 2 K X p 2 n 3 pax 

P2 = -r, — 5p 2 + -r. — 5pi H 5a (30) 

p 2 K 2 p 2 n 3 pa 2 

The other conjugates II and p have been related to pi and p 2 in (fTU|) and are therefore 
also fully determined 

5T = ^ r 5a+—50+—5p (31) 
C v pa p 

n = —5cf>+^—5p+—5a + w-v + w 2 (l-2cf>) (32) 

PK^ p 2 K w pCHcf, 

p = -J—Sp+—5^+—5a + w 2 4>(l-^) (33) 



p 2 n fl pn n pa p 



with 



a<f> 1 = «i 1 - a 2 1 (34) 

a~ l = ^ar 1 + (1 - (f))^- 1 (35) 

Kf 1 = /tf 1 + k 2 - 1 - 2k z ~ 1 (36) 

K^ 1 = (pK^ 1 - K 2 ~ l + K- A ~ l (37) 

s" 1 = Ar 1 - (i - 0) 2 k 2 _1 (38) 
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Eqs. (j2HHSni) as well as (|31H33|) contain 6 static susceptibilities as compared to 3 in a 
1-fluid description. In addition to the specific heat Cy there are 2 thermal expansion 
coefficients (since there are 2 densities) and 3 compressibilities (2 diagonal and one cross 
term). Eqs. (j28ll3T)|) are linear in the deviations from equilibrium, while ()31fJ3"3"|) explicitly 
contains nonlinear corrections involving velocities. Of course, the coefficients can still be 
phenomenological functions of the scalar variables (e.g. T or a, p or p, p\ and p 2 and even 
w 2 ) giving rise to additional nonlinearities that come with (usually) small coefficients. 
Note that neglecting cross-susceptibilities either in (j2HHS01) or in (J3THS3J) denotes two 
physically distinct (and incompatible) approximations, the justification of either one is 
not obvious a priori. 

Of course, there are situations where one has to go beyond the approximation used in 
the static equations above. Describing spinodal decomposition of the fluids, e.g. by an 
energy density e = a (f 2 + b [(1 — 4>) ln(l — <fi) + x0(l ~~ 0)] + C( P 4 + d (Vi</>) 2 immediately 
leads to nonlinear and gradient terms w.r.t. 0. 

What is left is the determination of hi = hi — Vj^ij in terms of cross-couplings to 
other variables are not possible due to symmetry. Thus this part of the statics is identical 
to that of ordinary nematics and can be taken over without any change 

tftf = K ijkl Vin k (39) 
hi = 6£-2^(V 3 ^)(Vw k )-Xa(H-n)H i -e a (E.n)E i (40) 

with Kijki = KxSjjSu + K 2 n p e P ijn q e qk i + K 3 njni5^,, the Frank gradient energy, and 
the transverse Kronecker symbol, 8^ = 5^ — riiUj. Orientation effects due to static 
external magnetic and electric effects enter through the diamagnetic (xa) an d dielectric 
(e ) anisotropy. For positive anisotropics the director is parallel to the external magnetic 
or electric field in equilibrium, which leads to a restoring torque outside equilibrium, e.g. 
to a (linearized) contribution to hi = XaH 2 5rii (with riibrii = 0). For negative anisotropies 
the director is perpendicular to the external fields and e.g. hi = \x a \(H ■ 5n)Hi. 

Since and hi are proportional to p 2 , so are the K n 's (and Xa, e a)- Again is 
linear in the deviations from equilibrium, but the inherent dependence of the material 
tensor on the direction n leads to nonlinearities in (|40|). 



5 Phenomenological Part of the Dynamics 

We now close our system of equations by setting up the connection between the currents 
and the thermodynamic conjugates (or rather their gradients usually called thermody- 
namic forces). For the irreversible part this is done by writing the entropy production in 
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terms of the forces 

2R = KijiViT^VjT) + Aj(Vin)(Vjn) + 2D i f ) (v<n)(v j r) + 7^ ^ 

+ Vij k i{VjVi)(Viv k ) +€ , ij m i rrij (41) 

including heat conduction, diffusion and thermodiffusion (/c^-, Djj, , respectively, all 
of the form = K±8^ + K||Tiin_j), director orientational viscosity 71, and viscosity related 
to gradients of the mean velocity Vijki- The latter has a Uijki = v k Uj symmetry and is 
of the uniaxial form [2H] characteristic for nematic systems. The last term in (j4"Tj) de- 
scribes the mutual friction between the two species as will become clear below. In (|4~D 
we have neglected viscosity-like contributions involving the relative velocity Vj-mj, since 
there is already dissipation due tomj. A more complete discussion of viscosity in a 2- fluid 
discussion is given in the Appendix. The dissipation function given above is bilinear in 
the forces, an approximation commonly called linear irreversible thermodynamics. Nev- 
ertheless it leads to nonlinearities due to (implicit and explicit) dependences of transport 
tensors on the variables. 

According to the dissipative parts of the phenomenological currents then follow 
from differentiating R 

.(„<*.) = _ {dR)/{dVtT) = _ KijV]T _ p(f)il _ (j))d (p VjU (42) 

Y m = ( dR )/( d - hi )^^ 5 ±- h . (43) 

a t f s) = -(dR)/(dVjVi) = -u m Vrv k (44) 
= {dR)/{dm t ) = $ j m j (45) 
.(i, dis ) = _ (9jR)/(9Vj n) = -p^V j n-p0(l-0)d J f ) V,T (46) 



where we have introduced the usual form of the diffusion (Z?y = pdij) as well as the 



thermo-diffusion tensor (D-p = p<j)(l — (j))d^). The ratios d* /d* (with * e {_ 



(T) 

and d* /k* are called the Soret and the Dufour coefficients, respectively (the latter being 
neglected usually in liquids). The viscosity term in (|4*4"|) has the same form as in a 1-fluid 
description. For a more general treatment of viscosity-like contributions cf. Appendix. 

In ad-hoc treatments of 2-fluid systems the mutual friction of the two species is in- 
troduced via an interaction force f 12 in the momentum equations for the single fluids, 
P1V1 = and P2V2 = —fu preserving total momentum. The force is related to the 
velocity difference, f 12 = — £ pi p2 w and is non-zero only if both fluids are present. This 
translates directly into w = —j^pw and can be compared to ([45 j). First, in a nematic 
environment the force f 12 is not necessarily parallel to w due to the possible anisotropy, 
rendering the £ to be a tensor ^ = (,±5^- + £||7ij7ij. Then comparison with ([45)1 gives 
£y = 0(1 — 0)£L, which shows that the ad-hoc choice for f 12 is the only possible one 
within linear irreversible thermodynamics. Of course, there is room for suitable nonlinear 
extensions (e.g. ~ f\ 2 or £ being a function of scalar state variables like T, 0, or p etc.). 

The reversible part of the dynamics is either dictated by symmetries or phenomenolog- 
ical. The symmetry parts have been discussed in secEland are listed in eqs. (|5IH24*j) . The 
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phenomenological reversible currents cannot be derived from any potential (especially not 
from any kind of Hamiltonian, despite being reversible). They are most easily derived by 
writing down all symmetry-allowed contributions to the various currents and then make 
sure that the entropy production (f27)J) is zero. We find 

Y^ rev) = -X ijk VjV k - X^VjTUk + PimjVjUi (47) 
a-- ev) = -\ kj i h k + 2(3 2 mi wj + fyimgj + mj9i) (48) 

xt v) = Pij v// + 7 « v,n _ Vj (\$ //,) • i. ^ V t n 3 

+(/3 2 wj + P 2 gj)(VjVi + ^i v i) + forrijiVjWi - ViWj) 
+p4Wj(VjVi - ViVj) + ^{niimjVj - m 2 Vi)F (49) 
jf' rev) = fam, (50) 

. (l,rev) /r 1 \ 

Ji = li0 m 3 ( 51 ) 



with 2X ijk = Xi8± n k + X 2 8^ rij, 2X^ = \[ m '8± n k + X 2 8^ nj, and /% = $±8+- + (3\\ 
and jij = 7j_5^ + 7|| n^Uj and F any function of the scalar variables or conjugates (e.g. T, 
p, 0). Since the term involving F is already of cubic order, we will neglect it in the follow- 
ing and suppress similar terms in other equations. Of the four flow alignment parameters 
A, only three are independent as will be discussed below. The /^-tensor in (|49|5()jl de- 
scribes a reversible entropy (energy) current due to a non-zero velocity difference as well 
as a change in the velocity difference due to a temperature gradient. In the limit of large 
£ the 7 and (3 parameters are related to diffusion and t her mo diffusion (see below). The 
1-fluid description is obtained in the limit of very large mutual friction, £ — > oo, which 
implies w — > 0. The usual nematodynamics (with an additional concentration variable) is 
regained, while p£w stays finite accommodating Eqs. (|19|45|49jl and is slaved by the other 
variables. 



6 Convective Velocities, Stress Division, and Con- 
centration Dynamics 

In (j47fJ5*T]) we have introduced terms, which are compatible with symmetries and R = 0, 
involving quadratic nonlinearities in the different velocities. Among them the f3\ term has 
a form quite similar to the convective term in (|20|) . Thus the actual velocity, with which 
^ is convected, is v conv = p~ 1 pi{l3\p 2 + + p~ 1 p 2 (—^iPi + 1)^2 an d can be either 
v±, v 2 or something in-between, depending on (3\. Since it is hard to imagine that rij is 
convected with a velocity larger than max{\ V\\, \ v 2 \), (3i is bounded p^ 1 < f3\ < —p 2 - 
If one accepts the reasonable assumption that n, is convected with the velocity of the 
nematic fluid v 2 , then (5\ = —p 2 l is fixed (while for (3\ = 0, there is v conv = v). 
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The choice of v conv has additional implications for the flow alignment parameters A. 
Since the director does not rotate in a frame that corotates with it, the quasi current 
couples to the vorticity by Y{ ot = eiji-njUjf; onv , where 2uj£° nv = curl-Uco^. For v conv = v 
this implies X^ = A 2 an d A2 — Ai = 2, or in the usual parameterization Ai = A — 1 and 
A2 = A + 1. For v conv = v 2 the conditions are 

Ai = A-l, A 2 = A + 1, -p 2 A< m) = A-l, -p 2 A 2 m) = A+l (52) 

In any case only two of the A's are independent. 

The phenomenological contribution ~ 7^ in (|HT|) affects the convection of the densities 
pi and P2 in (jEJUEJ). For 7_i_ = 7y = the densities are convected with v\ and v 2 , respec- 
tively, while (ft in (JT5|l moves with the velocity (l/p)(p 2 t>i + p±v 2 ). All these quantities 
G°ij Pi, Pi an d <p) are convected with v for 71 = 711 = — 1. Of course, other choices of the 
7's (made either by a theorist or by nature!) will lead to different convection velocities. 

Somewhat more involved is the question of the convective velocity for the different 
momenta. Because of g = pv, the total momentum has to be convected with v requiring 
/3' 2 = 0. If also w is convected with v then /3 2 — /?3 = Ai = 0, additionally, with the 
consequence that also the individual momenta, g l and g 2 , are convected with v. On the 
other hand, for g l and g 2 to be convected with 13 \ and v 2 , respectively, i.e. 

ft (1) +V j ft (1) ^ (1) +X« = (53) 

= (54) 



where 



X} 1] = PlElXi + ^Vjaij-^WiVjmj-mjVjvP (55) 

X/ 2 -* = — t-^-Xi + —VjOij — —WiVjirij + mjVjvf- 2 ^ (56) 
p p p 

in order to be compatible with ()18I19|) . the /3-parameters have to be (3 2 = 1/2, j3 3 = 
(1/pi) — (l/p 2 ), and /3 4 = 1/2 , thus ensuring that Xj and X 2 do not contain additional 
transport terms. This choice of parameters results in (P2U1 + p\V 2 ) / p to be the convective 
velocity for w (which is the same as for 0, when pi j2 are convected with i>i, 2 ). In addition 
the momentum current density due to flow then reads gjVi+Wjrrii = piV^Vj 1 ^ +p 2 v^Vj 2 \ 
which is the expected expression. 

The terms proportional to hi in (148149(1 constitute forces due to the nematic orienta- 
tional elasticity. Generally they act on both fluids. Using (|55|5fiJ) they read in linearized 
form 

9i l) \nem = <P (hji + P2^kji) ^ jh (57) 
9i \nem = (1 ~ 4>) {^kji ~ Pl^kji) ^ftk (58) 

Hence, for = (X^ = = X 2 m ^) this nematic force is distributed on fluid 1 and 
fluid 2 according to the ratio of pi/p 2 . It should be noted, however, that this kind of 
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nematic stress division is only compatible with the choice of n being convected by v, 
while it is incompatible with the choice of v 2 as the convective velocity for n (cf. (I52j) ) . 
Another reasonable case for the stress division problem is obtained for A^ = —P2^ijk 
(Ai = — P2A1 and A2 = — P2A2 ). I n ^at case the nematic force only acts on fluid 2 
(the nematic component). This case is compatible with fluid 2 being convected with v 2 
and the force then reads 

gP \nem = £ ((^ ~ ^fcj' 71 * + (A + 1)8^71 j)Vjh k = X k jN jh k (59) 

which is the form familiar from 1-fluid nematics. Thus, the so-called stress division 
problem (how hi in the stress tensor is divided between the two fluids) depends not only 
on specific material properties expressed by the phenomenological parameters and 
Ai )2 , but is also intrinsically linked to the question of the appropriate convection velocity. 

In order to verify experimentally the choices above, it seems to be difficult to directly 
measure specific convection velocities or the division of the nematic stress. However, there 
are situations, where these choices can be verified indirectly. Linearizing the dynamic 
equation for the relative velocity (|19J22J45J49j) and Fourier transform it w.r.t. time, w 
can be expressed by all the other variables. This can be used to eliminate w e.g. from 
the linearized dynamic equation for the concentration (JT5j) leading to 

iujcj) - r/;/ r V,V,H - ^4 T)e//V * V i T + AW (™ • ^) divh = (6°) 

where contributions of order 0(V 4 ) have been neglected. The effective diffusion and 
thermo-diffusion (Soret) coefficients have got additional frequency dependent contribu- 
tions due to the 2-fluid degree of freedom 

it „ = dr+ e^ii+j^ (61) 

p z pt,* + iu 

d (T)eff = ,(T) + &(l + 7*) (62) 

p + iu 

where the subscript * stands for either n or i. There is also a dynamic coupling to the 
nematic degree of freedom due to 

2p 2 V 1 ptw+zu 2 pU + iu) v ; 

These possible additions to the concentration dynamics, however, depend on the choices 
for the convection velocities as well as on the way how the nematic stress has been divided 
among the two fluids. Assuming the densities p\^ to be convected with the mean velocity 
v (implying 711 = 7j_ = —1) the additional contributions to diffusion and thermo-diffusion 
are all zero as well as the coupling to the nematic director. On the other hand, for p 1)2 to be 
convected with v\p respectively (7y = 7± = 0), both, diffusion and thermo-diffusion show 
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a dispersion step around lu ~ p£. For a nematic stress division among fluid 1 and 2 ac- 
cording to the ratio p\j p-i (implying X^ = A 2 = 0) there is no dynamic influence of the 
nematic degree of freedom on the concentration (A^ = 0), while for any other choice there 
is one. In particular, if only fluid 2 carries nematic stress (and p± t 2 are convected with Vx t 2, 
respectively), this dynamic coupling is given by A^ = — (pi/2p 3 ) [A^ 1 + £7/) — C^ + ^I 1 ] 
for strong friction (p£* ^> uj). 



7 Simplified 2-Fluid Nematic Equations 

In the preceding sections we have derived the most general and complete set of 2-fluid 
equations for a nematic and Newtonian mixture. Special emphasis has been laid on the 
correct form of the nonlinearities that come with the 2-fluid description. However, these 
equations are for most purposes unnecessarily complicated and can be simplified using 
reasonable assumptions. Starting from the correct general equations such assumptions, 
clearly spelled out, lead to controlled approximations and to a set of 2-fluid equations, 
whose limitations and implicit assumptions are clear and well defined in contrast to most 
ad-hoc approaches. 

Here we want to display explicitly 2-fluid hydrodynamics for a nematic/simple fluid 
mixture under the following assumptions, 

a) convection with natural velocities (for n, g 2 , p 2 and g 1 , p\ this is v 2 and v±, respec- 
tively, or explicitly (5\ = — p 2 x , P2 = Pa = 1/2, /3 3 = pT/ 1 — p 2 A ", J± = = 7y and eq. (|52*j) ): 

b) the linearized orientation-elastic force acts on the nematic fluid (index 2) only (i.e. 
Ai,2 = — P2A1™ ); 

c) global incompressibility, 5p = (i.e. Spi = —Sp 2 ); 

d) neglecting the phenomenological reactive entropy current (/3jj = 0); 

e) linearizing the phenomenological dissipative currents, but keeping quadratic nonlinear- 
ities otherwise. 

Then the following set of equations is obtained: 
The incompressibility condition 

= divv (64) 
or = w ■ 'V p 1 + pidivvi + p 2 divv 2 (65) 
or = w • V0 + <p div(l — <p)w — (1 — <p) divcpw (66) 

the concentration dynamics 

+ Vi ((jwi + 0(1 - <f>)wi) - dijViVjipi - p 2 ) ~ 0(1 - fldjpVjViT = (67) 

or pi + vi- Vp 1 + p 1 divv 1 -pd lJ V i y j (p 1 -p2)-—d(pV i V J T = (68) 

p 

or p 2 + v 2 - Vp 2 + p 2 dwv 2 + pdijViVjip!- p 2 ) + —dpViVjT = (69) 
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the entropy dynamics (heat conduction equation) 



the nematic director dynamics 



hi + v^Vjrii - \ijkVjV^ 2) - —Xi jk WkVj4> + 7i hj = 0. (71) 

P2 

There is a (nonlinear) coupling to the concentration variable, which is not possible in the 
1-fluid description. 

For the momentum balance of the two different species we get 

ptvV + pxvPVjvV + ^i Vi (p + \p 2 {v\ -v 2 2 )) + ^Vi(^i - /2 2 ) + -V^Vm) 

p 2 ' p p 

+ —hjVirij + —X kj ih k Vj4> + iijPiP2Wj - K^kWj^i 11 ^ = ( 72 ) 

P P2 

P2VP + PlV^VjvV + ^V,(p - \ Pl (v\ - V\)) - ^V,(^ - p 2 ) + ^V^ViUk) 

p 2 - p p 

- —hjVirij - —\ kji h k Vj<f) - Vj(X kji h k ) - iijpxpiWj - vSuVjVivj® = (73) 

P P2 

Note that although we made the approximation that the linear orientational-elastic stress 
does only act on fluid 2, there are inevitably nonlinear contributions to fluid 1, too. There 
is also a (nonlinear) coupling of fluid 1 to the concentration, if nematic distortions (hi ^ 0) 
are present. In ()72J73j) cross- viscosities have been neglected (cf. Appendix). 

In order to facilitate actual calculations we also give eqs. (|72|73[) as dynamic equations 
for the total momentum and for the relative velocity 



pit + Vip + pVj (viVj + 0(1 - <j))wiWj \ + Vj (^/ kj Vin k - \ kji h 



i'l'k 

VijkNjViVk = (74) 



Wi + (vj + (1 - 2(j>)wjJ VjWi + Vj (pi - fjQ + v • w + - (p)w 2 j + p&jWj 

+—h J V l n J +V J (—X kji h k ) = (75) 

P2 V P2 ' 

Note that the approximation for the viscosities made in (|74I75|) is not compatible with 
that used in ()72I73|) ; their interrelation is discussed in the Appendix. 

Due to the incompressibility condition the pressure is no longer an independent vari- 
able nor is it given by the other variables (i.e. (f5f or cannot be used), but it serves 
as an auxiliary quantity to ensure the incompressibility for all times, i.e. divi; = 0, which 
leads to the very complicated condition that determines Sp 

Ap = -V % V 3 { Pl vf ] vf ] + P2 vf ) vf ) )-V l V j ^kF l n k ) (76) 
+AVjVj (5 kj nih k ) + VijkiViVjViv k 
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Although Sp does not show up in the dynamical equations, it is still present in boundary 
conditions etc. and it contains combinations of the viscosities different from those present 
in the incompressible dynamical equations. 4 In contrast to 1-fluid descriptions for simple 
fluids, where the incompressibility condition leads to a considerable mathematical simpli- 
fication, this is no longer the case for a 2-fluid description due the complicated form of 
(|76|). even if incompressibility is a very good approximation in physical terms. 
Of the statics (J2EHHH1) only the following equations remain 

ST = TCy 1 5a + a^5<f> (77) 

with 5(p = p^Spi = —p~ 1 5p 2 , while (JHHJ) and pHjl are unchanged. Note that 5fi is not 
needed, but follows from 5p via eq. ©• 



8 Isotropic Viscoelastic Fluids 

In this section we discuss the 2-fluid description of isotropic viscoelastic fluids by choosing 
a Newtonian fluid as fluid 1 and an elastic medium as fluid 2. The latter can be a 
permanent network (showing e.g. diffusion) or a temporary one relaxing on a finite time 
scale. The considerations for setting up a complete nonlinear 2-fluid description for such 
systems is quite similar to that of the 2-fluid nematic discussed in detail in the previous 
sections - only that the nematic degree of freedom VjUi has to be replaced by the Eulerian 
strain tensor Uij, which we use to describe the elastic degree of freedom. In the following 
we present an abbreviated discussion, starting with the general energy expression. 

8.1 Thermodynamics 

In analogy with the development in section^ the general energy expression for an isotropic 
elastomer network immersed in a Newtonian solvent is given by 

E = tV = j edV = E{M U M 2 , V, G 2 , S, M 2 U ij ) (79) 

from which the conjugate quantities are derived. The elastic stress, conjugate to the 
strain, is $^ = dE/d(M 2 Uij) = de/d(p 2 Uij) = p^&ij, while the definitions of the other 
conjugates (j2J) remain unchanged (except that they are to be taken at constant rather 
than constant Vj-rij). With the new definition /2 2 = fJ>2 + p~ x ^ijUij the expressions for the 

4 Incompressibility 'reduces' the number of independent components of the viscosity tensor from five 
to three (in the uniaxial case), only if a redefinition of the pressure is done, cf. JJJJ] p.41f and [33] 
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pressure (|3I7|) and the relations of the different sets of conjugates (fTU|) remain unchanged, 
while the Gibbs and Gibbs-Duhem relations read 



de = Tda + pi dpi + p 2 dp 2 + v 1 ■ dg 1 + v 2 ■ dg 2 + ^ dUij (80) 

= T da + II d(f) + p dp + v ■ dg + m ■ dw + $^ dUij (81) 

dp = a dT + pi dpi + p 2 dp 2 + g 1 ■ dvi + g 2 ■ dv 2 — dUij (82) 

= a dT + p dp + g ■ dv — m ■ dw — $y dUy (83) 

Rotational invariance of the Gibbs relation ()80I81|) leads to the condition 

Uik^kj = U jk <S> ki (84) 
which, as is seen later on, ensures the stress tensor to be symmetric. 

8.2 Statics 

The conjugate quantities defined by the Gibbs relation (j8QI81|) are linked to the variables 
by a set of phenomenological equations containing static susceptibilities as parameters. 
This constitutes the static part of the hydrodynamics. Instead of the nematic molecular 
fields hi and ip^ we now have the elastic stress $jj as conjugate field. As a symmetric 
2-rank tensor it consists of a scalar quantity, the trace and the deviator = 
$ij — (l/3)5ij&kk- Being a scalar & kk can couple to the other scalar variables like densities, 
concentration or entropy, just like the 3 other scalar conjugates {T, II, p or pi, p 2 } by 

5T = ^5a + —5 Pl + —5p 2 + —U kk (85) 

Cy pOL\ pa 2 Ct3 

T r 1 SJ . 1 r 1 rr 

= d0H d)pH U kk 

C v ol^ pa p a 3 

1111 
pi = — — 5p 1 + — — 5p 2 H U kk (86) 

p K\ P &3 POCI pK,^ 

II = — 5<f) +— 5p+ — 5a + — U kk + w-g + pw 2 (l-2<f)) (87) 
P2 = -^—$P2 + -^—$Pi + —5a + —U kk (88) 

P K>2 P K3 P«2 P«5 

1111 

p = ~^5p + 5(f) + 5a + U kk + w 2 (p(l - 0) (89) 

P Z K /1 pK n pCXp pK p 

111 

<$>kk = ciU kk + —5a + — 5p! + —5p 2 (90) 
= ciU kk + 

c 

= cUU^-^Ukk) (91) 



«3 


pn± 


pn 5 


1 - 




1 . 


— 5a -\ 


- 1 5(f) 


+ 5p 


a 3 


P^u 


pn p 
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where - in addition to (jMHHEI) 



Kp~ l = <\)Kf x + (1 (92) 
k u ~ 1 = nf 1 - K 5 _i (93) 

involving 2 new generalized compressibilities K4 5 or k U)P and one expansion coefficient 
related to the trace of the elastic strain U kk . It should be noted that for real solids at 
finite temperatures U kk 7^ <W P i n contrast to ideal elasticity theory. The reason are the 
point defects, which allow not only the dissipative motion described above, but also static 
temperature and pressure changes due to Ukk even at constant density The new static 
susceptibilities q and c tr are the usual elastic moduli of Hooke's law. The longitudinal 
one is related (in addition to the compressibility k m ) to the sound velocity. The transverse 
modulus leads to transverse sound, which is however relaxing due to (jl04|) if Q and Q r 
are not zero. 



8.3 Dynamics 

The dynamical equations for the elastomeric and solvent degrees of freedom are 





= 


(94) 


<j> + VjVrf + p~ l Vi (p0(l - <f>)wi + jty 


= 


(95) 


e + Vj{e + p)vj + Vij} e) 


= 


(96) 


■ 1 T7 / 1 ■ (o~,rev) ■ <a,dis) \ 
<T + Vj{<TVj+3> '+]> ') 


= R/T 


(97) 


wi + vjVjWi + v,n + x} rev) + x} dls) 


= 


(98) 


k + VjQiVj + Vip + Vji-Qij + $ jk U ik + <$> ik U jk + a£ ev) + a}f s) ) 


= 


(99) 


Uij + v k V k Uij + U kj ViV k + U ki VjV k - + Z^' ev) + z}f ls) 


= 


(100) 



the first 5 equations have the same form as before (but different phenomenological cur- 
rents, see below). In the dynamic equation for the strain ()100|) there are nonlinear cou- 
plings to the velocity gradient that have the form of the so-called lower convected deriva- 
tive j3H ESI an d the appropriate counter terms show up in the stress tensor as additions 
to the phenomenological parts. In all dynamic equations the convective velocity chosen is 
the mean velocity v, since this allows a simple and thermodynamically consistent way of 
writing the equations. However, as in the case discussed previously there are phenomeno- 
logical terms in the reversible currents that allow a different choice of the convective 
velocities (see below). 
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8.4 Currents 



Following the previous development, we can establish the reversible and dissipative cur- 
rents, and their constraints. 

For the phenomenological parts of the currents there is the condition 

R = -j^ViT + IT V^f - aVVjVi + QijZV + m, X« > (101) 

with the equal sign (> sign) for * = rev (* = dis), respectively. 

The dissipative parts of the currents introduced above can again be deduced from a 
dissipation function that reads in bilinear approximation 

2R = K(VT) 2 + D(Vn) 2 + 2D (T) (VT)-(Vn)+£'m 2 

+(ijkl®ij®kl + iijklmn{^i^jk){^l^mn) + V ijk l{V >i)( V 

((V^)(V,m fe ) + (V J m J )(V^ fc )) + ^(V^V,™*) (102) 

where all 4-rank material tensors have the form = Ui5ij5 k i + (l/2)v tr (5 ik 5ji + 5u5j k — 
(2/3)6ij5ki)) and (,ijki mn contains 4 parameters £1-4. In the C-tensor Q and ( tr are describ- 
ing the relaxation of elastic strains and the £1-4 give rise to vacancy diffusion as can be 
seen in the following expressions 

jt 4is) = -(dR)/(dV i T) = -KV i T-p<f>(l-<f ) )d^V i n (103) 

Z^ s) = (dR)/(d$ ij )=t m Qki-Vk(Zw m nV& mn ) (104) 

a (dis) = -(dR)/{dV j v i ) = -v m V l v k -v$ l V l m k (105) 

x m = { d R)/{dm l )=i'm l -V {v^ l V l w k + ^ kl V l v k ) (106) 

j(,i,<us) = -( 5j R)/(9v l n) = -prfv J n-p0(i-0)rf( T )v J T (107) 

where diffusion and thermodiffusion is written in the usual way with D = pd and = 
p4>(l — (f))d^ T \ For a permanent network that does not relax, the relaxation parameters 
vanish (Q = ( t -r = 0). For the reversible parts of the currents we find 

z}- ev) = \ {u) {ViTrij + VjUii) + f3 7 {U k] V l m k + UkiVjmk) + /3 6 m fc V fc £^ (108) 
a% ev) = 2(3 2 m iWj (109) 

X} rev) = 2V j (AW$ ij ) + WiT + 7 V,n - p 6 $ kj ViU kj + p 2 WjiVjVi + V iVj ) 

+V j f3 7 ($ kj U ik + <& ki U jk ) + PsmjiVjWi - V^) + /^(V^ - v>j) (110) 

jt' rev) = firm (111) 

j?' rev) = im (112) 

8.5 Convection, Stress, and Concentration Dynamics 

As in the case of the 2-fluid nematics the velocities which with the variables are convected 
can be tuned by choosing special values for the coefficients P n and 7. E.g. for 7 = 



19 



the densities pi j2 are convected with respectively (and the total density p and the 
concentration <ft with v and (1/ p){p 2 Vi + Pi v 2)-, respectively), while for 7 = — 1 all 4 
quantities are convected with v. Similarly, for /5 4 = 1/2 = (3 2 and j3 3 = (1/pi) — (I/P2) 
the momenta gi ; 2 are convected with vi i2 (and the total momentum g and the relative 
velocity w with v and {1/ p)(p 2 v\ + P1V2), respectively). For /3 6 = — l/p 2 the strain L^- 
is convected with v 2 and for f3 7 = — l/p 2 the lower convected derivative contributions 
in (I100|) effectively come with u 2 (producing an additional cubic term in ()110|) ~ Vjp 2 , 
which can be neglected as other cubic terms). Even the convection of the entropy can be 
tuned by choosing (3 = /3 + j3 Q0 a where /3oo = 1/pi, = 0, = — l/p 2 leads to the convective 
velocity to be Vi, v, v 2 , respectively. 

The distribution of the elastic stress among the two fluids is governed by the coefficient 
A^. For, respectively, 2\( u > = l/p 2 , = — 1/pi, or = 0, the elastic stress is carried by fluid 
2, fluid 1, or is equally distributed between them. 

As in the case of 2-fluid nematics we can linearize and Fourier transform the dynamic 
equations, thus eliminating w from e.g. the concentration dynamics. Neglecting fourth 
order gradient terms we get 



iu<f) - d eff AU - ^d {T)eff AT - 2A w V i V i $y = (113) 



with frequency dependent effective diffusion and thermo-diffusion coefficients 



= rf+^^ T (114) 



P1P2 (7 + I) 2 
p 2 p^ + iuj 

d (T)eff = d (T) + gT + l) (n5) 

p^ + iuj 



and the dynamic coupling to the elastic degree of freedom by 

A W = PiPl X (U)l+JL (116) 
p pi + iuJ 

Again these possible additions to the concentration dynamics, however, depend on the 
choices for the convection velocities as well as on the way how the nematic stress has 
been divided among the two fluids. Assuming the densities pi )2 to be convected with 
the mean velocity v (implying 7 = — 1) the additional contributions to diffusion and 
thermo-diffusion are all zero as well as the coupling to the nematic director. On the other 
hand, for pi j2 to be convected with t>i i2 respectively (e.g. 7 = 0), both, diffusion and 
thermo-diffusion show a dispersion step around u ~ p£. For the elastic stress division 
among fluid 1 and 2 according to the ratio p\j p 2 (implying A^ = 0) there is no dynamic 
influence of the elastic degree of freedom on the concentration, while for any other choice 
there is one. In particular, if only fluid 2 carries elastic stress (A^ = I/P2), this dynamic 
coupling is given by A^ = pi/p 2 ^ for strong friction (p£ ^> iS). 
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8.6 Simplified elastomeric two-fluid equations 

In the preceding sections we have derived the most general and complete set of 2-fluid 
equations. These equations are for most purposes unnecessarily complicated and can 
be simplified using reasonable assumptions. Starting from the correct general equations 
such assumptions, clearly spelled out, lead to controlled approximations and to a set of 
2-fluid equations, whose limitations and implicit assumptions are clear and well defined 
in contrast to most ad-hoc approaches. 

Here we want to display explicitly 2-fluid hydrodynamics under the following assump- 
tions, 

a) convection with natural velocities (for C/^, g 2 , P2 and g l: p 1 this is v 2 and v±, respec- 
tively, or explicitly (3 7 = -p 2 l = A>, fa = & = 1/2, f3 3 = p\ l - p 2 l , 7 = 0); 

b) the linearized elastic force acts on the elastomeric fluid (index 2) only (i.e. 2A^ = p 2 r ); 

c) global incompressibility, bp = (i.e. 5p\ = —Sp 2 ); 

d) linearizing the phenomenological dissipative currents, but keeping quadratic nonlinear- 
ities otherwise. 

Then the following set of equations is obtained: 
The incompressibility condition 

= dlvv (117) 
or = w ■ Vpi + pidiv V\ + p 2 div v 2 (118) 
or = w ■ V0 + 0div(l - 4>)w - (1 - 0) div0it> (119) 

the concentration dynamics 

+ V, (<M + 0(1 - 0H) - c^ViV^pi - p 2 ) - 0(1 - 0)4 T) VjV,T = (120) 
or p 1 + v 1 -Vp 1 + p 1 divv 1 -pd ij W i W j (p 1 -j2 2 )-^dj ; pV i W j T = (121) 

or p 2 + v 2 -Vp 2 + p 2 divv 2 + pd i ,V i V j (p 1 -p 2 ) + ^d^ ) V l V j T = (122) 
the entropy dynamics (heat conduction equation) 

a + v i V l a + ^VifafrWi) - /^V^T - —d^V \V - p 2 ) = (123) 
the elasticity dynamics 



p p « 



Uv + v^VMj - -(VjvP + Vivf } ) - |(^V, + W&) In ^ + U ki V ^ + U kj V^ k 2) 

+ CiSij&kk + (tr{<5>ij - ^SijQkk) - ZiSijAQkk - - ^(ViVj^fcjfc + ^V fe V^ w ) 

- ^(ViVfc^-fc + VjVA) = (124) 

There are nonlinear couplings to the concentration variable (the cubic one has been sup- 
pressed), which are not possible in a 1-fluid description. 
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For the momentum balance of the two different species we get 

Piv^ + PivWVjV® + ^Vi(p + \p 2 {v\ - vl)) + ^V<(//i - Jh) + -$kjViU kj 

P I ' p p 

- p&yVi In ^ + P -^(3V{T + Z ijPl p 2Wj - vSjfrjVtvj? - vSSVfrv® = (125) 

p p j j 

•(2) i (2)v7 ( 2 ) i P 2 V7 I 1 I 2 2\\ PlP? V7 C - ^ ^ * Y7 7T 

' p 2 p p 

+ p&jVj In ^ - ^(3ViT - + VjiQjkU* + $> ik U jk ) - ^pip 2 ^ 

- 4W;V^ fe (2) - vSSVjVtvP = (126) 

Note that although we made the approximation that the linear elastic stress does only 
act on fluid 2, there are inevitably nonlinear contributions to fluid 1, too. There is also a 
(nonlinear) coupling of fluid 1 to the concentration, if elastic distortions are present. 

The different approximations for the viscosities are discussed in the Appendix. 

In order to facilitate actual calculations we also give eqs. (|125ll26j) as dynamic equa- 
tions for the total momentum and for the relative velocity 



pvi + ViP + Vj (pViVj + ^^-WiW^ - 



PXP2 (c) 

p u m 

Wi + (vj + P2 - Pl W^jVjWi + Vj (p\ - p 2 + V ■ W + — ' ^ 1 W 2 ^j + p£ijWj + Vj — Q>ij 

+ -^ kj ViU kj - -Vj^Ua) - ^v$V&jw k - v^ j V l v k = (128) 

p2 P2 P 

In order to conserve the global incompressibility condition for all times, i.e. divi; = 0, 
the pressure has to fulfill the relation 

Ap = -V t Vj (pM 1} vf + P2 vf ) vf ) ) + ViVfin - ViVi {® kj U ik + <S> ik U jk ) 

+v ijk N l V 3 V l v k + p l p 2 p- 1 v$,y i V j V l w k (129) 

In contrast to 1-fluid descriptions for simple fluids, where the incompressibility condition 
leads to a considerable mathematical simplification, this is no longer the case for a 2-fluid 
description due the complicated form of ()129J) . even if incompressibility is a very good 
approximation in physical terms. In particular, Ap is not only connected to compressions 
(U kk ), but also to shear deformations, even in linear order. 
Of the statics (jHSHSDI) only the following equations remain 

ST = TCy l 8a + a^5(t) + a^ l U kk (130) 
8(pi-p 2 ) = P" 1 ^ 1 $4> + P -1 "^ 1 £cr + K7 a l U kk (131) 
§ kk = ciU kk + a^Sa + p' l Kl l 5<p (132) 



+ 2V j {$ jk U ik ) - v m VjVtv k - —v&VjViWk) = (127) 

p 
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with 5(fi = p 1 5pi = — p 1 Sp 2 , while eq. ()9ip remains unchanged. Note that 5/j, is not 
needed, but follows from 5p via eq. (JH3) • 



9 Discussion 

Within the general framework of hydrodynamics and thermodynamics we have set up 
a consistent nonlinear 2-fluid description of complex fluids, in particular for lyotropic 
nematic liquid crystals and polymer solutions or swollen elastomers. Such a general 
theory determines the frame for any ad-hoc model, which has to be a special case of the 
general one. The comparison with the general theory also reveals implicit and explicit 
assumptions, approximations and possible generalizations of a given model. A simple or 
" natural" choice in a given model may not be mandatory, but rather imply a presumption. 

Quite generally we find that neither the velocity, with which a certain variable is 
convected, nor the stress division between the different fluids can be determined by general 
principles, but is rather system or material dependent. On the other hand, there are 
certain restrictions and interrelations among the convective velocities and other physical 
effects that limit the possible choices. For the two densities pi, p 2 e.g., the natural choice 
for the convection velocities seems to be their native velocities V\ and v 2 , respectively. 
This implies that the total density is convected with the mean velocity v (as required by 
mass transport), while the concentration is convected with (1/ p)(p 2 V\ + P1V2). Another 
obvious choice would be the mean velocity as convection velocity for both, the total 
density as well as the concentration implying that also p% and p 2 are convected with v. 
However, the actual convection velocity depends on the value of the material dependent 
(reactive) flow parameters 7j_ and 7y, defined in eq. (j5Tf . 

For the nematic degree of freedom the convective velocity again depends on a material 
parameter (/3i defined in eq.flUJ)) an d is not necessarily equal to v 2 (if fluid 2 is the 
nematogen). However, the value of (3\ influences also the flow alignment of the director 
(and the back flow due to director reorientation), which can be measured in shear flow 
experiments. In the case of visco-elastic and elastic media, which are described by a 
dynamic equation for the (Eulerian) strain tensor Uij, there are two velocities involved. 
One is the usual convection velocity (vkVkUij) and the other one occurs in the "lower 
convected" part {U^jViVk + UkiVjVk). There is no fundamental reason for the two to 
be equal and their actual value depends on the (reactive) flow parameters f3§ and f3y, 
respectively, defined in eqs. (|1()8[ IllOj) . 

For the evolution equations of the momenta special care has to be taken to get a 
description, which is compatible with general laws (cf. Chapter IHJ). The currents and 
quasi-currents that enter the description in terms of either the total momentum and the 
velocity difference or the two individual momenta are not the same as seen in eqs. fl55|l56|) . 
In the nematic case the stress division problem depends on the flow alignment parameters 
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as well as on the convection velocity of the director, while in the visco-elastic case the 
crucial material parameter A^, eq. (jl08lll0|) . is not related to a convective velocity. The 
delicate question of viscosities, and approximations related to them, is discussed in detail 
in the appendix. 

A prominent feature of the 2-fluid description is the coupling of the concentration 
dynamics to the velocity difference. This leads to a frequency dependent effective diffu- 
sion and thermo-diffusion, as well as a frequency dependent coupling to the nematic or 
the visco-elastic degree of freedom. For low frequencies these contributions to the con- 
centration dynamics constitute additional dissipation channels, while for the short-time 
dynamics (below the relaxation time of the velocity difference) they are reactive. 

Recently, 2-fluid descriptions of diffusion in polymeric systems have been given [HHl EZ1 
based on the GENERIC approach making use of Poisson brackets. A detailed comparison 
with these formulations is beyond the scope of this manuscript and will be discussed 
elsewhere. 



Acknowledgments 

This research was supported in part by the National Science Foundation under Grant No. 
PHY99-07949. 



Appendix 

In this Appendix we discuss viscosity and viscosity-like phenomena in the 2-fluid hydro- 
dynamics. We show that in order to get the 1-fTuid limit (the binary liquid) correctly, 
some care has to be taken when the usual approximations are made. 

If there is only one velocity present, the viscous contribution to the dissipation function 
is Vijki^J jVi)(y iv k ) with u ijkl = v jiU = u ijtk = u kUj , which ensures that only symmetric 
velocity gradients contribute to dissipation. Antisymmetric velocity gradients, curl?;, 
describe rotations. A solid body rotation (curl?; = const.), however must not increase 
the entropy and (curl?;) 2 contributions are not allowed in the dissipation function. With 
these symmetries the viscosity tensor has 2 coefficients for the isotropic 

2 

"ijki = "(Sjrfik + Sudjk - -6ijS k i) + C^ijSki (A.l) 
and 5 for the nematic case 

Vijki = v i {SjiSik + 5u5 jk ) + 2(ui + u 2 - 2u 3 ) ninjn k ni + (u 5 - v A + z^X^fc^z + 4z™i?v) 
+ (u A - u 2 ) 5ij5 k i + (zaj - ^{njUiSik + rij-nAi + nin k 5ji + UiUiS jk ) (A.2) 
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In a 2-fluid description the same restrictions hold with respect to the mean velocity v, 
since it is the conjugate to the momentum density and curl v = const, still describes solid 
body rotations. There are no such restrictions to the relative velocity w, and curlm can 
contribute to the dissipation. The most general form for viscous dissipation in a 2-fluid 
description thus reads 

2R^ = uwCVjViXyn) + 2i/$(V i « t )(V l m Jfc ) + ^ l (V j m i )(V l m k ) 

+ is}p (curl m)i (curl m)j + 2to^(curlra);V/z;& + 2z/A e fc ) (curlm) i V 3 -m fc (A.3) 

Note that only Uijki has the dimension of a viscosity, while v\f ki and vf^ k are kinematic 
viscosities, while i^), vfX, and v-p are viscosities divided by p 2 . The tensors v and 
z/W have the familiar form (|A.1J) or (|A.2|) . For there is no a-priori reason for a 
= v kiij symmetry, since v and m are not equivalent. However, as will be seen below, 

(c) 

a consistent 2-fluid description is only possible, if this symmetry holds and i/A^ has the 
form (|A.1J) or ()A.2j) . 5 The tensor ufp = u^Sij or v-p = V\5ij + u^TiiUj contains 1 
or 2 coefficients in the isotropic and nematic case, respectively. The 3rd rank material 
tensors, symmetric in the last two indices = v^f 1 are zero in the isotropic case and 
both carry one coefficient in the nematic case v^ k ^ = v^ d ' e \ei k injrii + €j k iriini). For the 
dissipative currents this leads to 

<? t f s) = -vw^iVk -Vijli^irrik - i/gj(curlm) fc (A.4) 
X- dls) = Cj m j - Vj (v^ k ) Vim k + 4t V ^ + ^ffefcitCcurl m)j + v$ k eijiV p v k 

+ v tl e iji^P m k + fjy*(curl m) fc J (A.5) 

Since there is already friction ~ m^, very often the viscosity-like dissipation ~ Vjirii is 
neglected altogether {y^ c > = = = u^ w > = = 0). Such an approximation leads 
to (|43I44|) . In the strong coupling limit, where w (and m) vanish, this approximation 
seems to be appropriate and it correctly gives the 1-ffuid limit of binary mixtures. On 
the other hand, for two fluids only gently coupled there is no a-priori reason, why e.g. 
the tensor v^ c > (or should be neglected compared to is, since both terms contain 

gradients of V\ as well as of v 2 - Indeed, the dissipation function in terms of V\ 2 reads 

+ ^t^ r) (curl [ Vl - ^ 2 ]) l (curl [ Vl - V 2 ]) j 

+ 2^(curl [ Vl - v 2 \H^V 3 vV + „<g> Vi ,f ) (A.6) 



5 Without this symmetry the form (jA.lfl still applies for the isotropic case, while in the nematic case 
an additional coefficient is present, i.e. the term ^(riinjSki + nkniSij) in (| A. 2|> splits into two different 
parts, u 5a mnj5ki + v 5b n k ni6 i3 . 
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where curl v is absent in the dissipation function. Comparing (jA.3|) and (jA.6j) in harmonic 
approximation, i.e. neglecting cubic and quartic terms involving e.g. (curlvi)^ Vj/>i, 
( Vl x VpO^V.px, or ^(ViPiXvrf^ ^^(^(V^), we get 

= piU W +pip2U ^ (A.7) 

= p 2 ^) - p lP2 v& (A.8) 

A$, = pl^ + 2 P 2 iP^% + P 2 ipHjii (A.9) 

= ^H-2pgpn/g, + p?^§2 (A.IO) 

2 (12) . / (c) (c) \ 2 2 M /a i i\ 

P "iW = PlWW + PlP2{p2Klij ~ Pl"ijkl) ~ PlP2 V ijkl (A.H) 

Again there is no a-priori reason for is^j = v^j~. However, since the tensors v, v^"\ 
and do have this symmetry, eqs. (jA.9|A.10j) force u^ c ' to have it, and finally (jA.llj) 
requires also z^ 12 -* to have this symmetry and thus to be of the form (jA.lJA.2j) . 

Neglecting the curl-terms means the same in both descriptions fjA.3j) and (jA.6j) . i.e. 
z/(' r ) = z/( a! ) = = v( d2 ) = 0. For the symmetric velocity gradient terms the ap- 

proximation = = used in (j44!45j) (i.e. no Vim^-terms in (jA.3jl ) leads to 
p|i/( 1 ) = p\v^> = \piP2^ 12 ^ leaving only one viscous tensor independent. Neglecting only 
the cross-viscosity in (jA.3j) does not imply the cross-viscosity in (jA.6j) to vanish, since 
z/( c ) = o gives pz/ 12 ) = p 2 v^ + pi^ 2 \ The opposite case z/ 12 ) = used in (j72!73j) leads to 
a non-zero p\piv^ = P2^^ — P\v^ (and v = + = p^[ 2 u^ + p~^ %v \ Thus, 

the approximations leading to (j44!45j) are not compatible to those used in (j72!73j) . 

In the general case (jA.6j) leads to the following viscous contributions to the left hand 
sides of (j72j) and J73J), respectively 

- ^^V,4" - "SSv^f - ^,g'v j( curl [„ - ,,1), 



and 



i/ 2) V V,?; (2) ;y (12) V V,?; (1) Plp2 iM^XJ (purl [7, 11 ]), 
u ijkl v j v l v k u klij V j v l v k p~ ki i V il curl ri v 2\)l 



'e^V^curl [ Wl - v 2 )) t + ^%,V, (i/g^V^ + ^Vpt^) (A.13) 



In the 1-fluid limit z/ dl > d2 ) have to vanish and /^z/^+z^ 12 )) — > pi^ and p(u^ +u ( - 12 ^) — ► p 2 z/, 
which is obtained for vanishing i/*") and z/( c ). 
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